117 research outputs found
Curl Constraint-Preserving Reconstruction and the Guidance it Gives for Mimetic Scheme Design
Several important PDE systems, like magnetohydrodynamics and computational electrodynamics, are known to support involutions where the divergence of a vector field evolves in divergence-free or divergence constraint-preserving fashion. Recently, new classes of PDE systems have emerged for hyperelasticity, compressible multiphase flows, so-called first-order reductions of the Einstein field equations, or a novel first-order hyperbolic reformulation of Schrödinger’s equation, to name a few, where the involution in the PDE supports curl-free or curl constraint-preserving evolution of a vector field. We study the problem of curl constraint-preserving reconstruction as it pertains to the design of mimetic finite volume (FV) WENO-like schemes for PDEs that support a curl-preserving involution. (Some insights into discontinuous Galerkin (DG) schemes are also drawn, though that is not the prime focus of this paper.) This is done for two- and three-dimensional structured mesh problems where we deliver closed form expressions for the reconstruction. The importance of multidimensional Riemann solvers in facilitating the design of such schemes is also documented. In two dimensions, a von Neumann analysis of structure-preserving WENO-like schemes that mimetically satisfy the curl constraints, is also presented. It shows the tremendous value of higher order WENO-like schemes in minimizing dissipation and dispersion for this class of problems. Numerical results are also presented to show that the edge-centered curl-preserving (ECCP) schemes meet their design accuracy. This paper is the first paper that invents non-linearly hybridized curl-preserving reconstruction and integrates it with higher order Godunov philosophy. By its very design, this paper is, therefore, intended to be forward-looking and to set the stage for future work on curl involution-constrained PDEs
A Hybrid Godunov Method for Radiation Hydrodynamics
From a mathematical perspective, radiation hydrodynamics can be thought of as
a system of hyperbolic balance laws with dual multiscale behavior (multiscale
behavior associated with the hyperbolic wave speeds as well as multiscale
behavior associated with source term relaxation). With this outlook in mind,
this paper presents a hybrid Godunov method for one-dimensional radiation
hydrodynamics that is uniformly well behaved from the photon free streaming
(hyperbolic) limit through the weak equilibrium diffusion (parabolic) limit and
to the strong equilibrium diffusion (hyperbolic) limit. Moreover, one finds
that the technique preserves certain asymptotic limits. The method incorporates
a backward Euler upwinding scheme for the radiation energy density and flux as
well as a modified Godunov scheme for the material density, momentum density,
and energy density. The backward Euler upwinding scheme is first-order accurate
and uses an implicit HLLE flux function to temporally advance the radiation
components according to the material flow scale. The modified Godunov scheme is
second-order accurate and directly couples stiff source term effects to the
hyperbolic structure of the system of balance laws. This Godunov technique is
composed of a predictor step that is based on Duhamel's principle and a
corrector step that is based on Picard iteration. The Godunov scheme is
explicit on the material flow scale but is unsplit and fully couples matter and
radiation without invoking a diffusion-type approximation for radiation
hydrodynamics. This technique derives from earlier work by Miniati & Colella
2007. Numerical tests demonstrate that the method is stable, robust, and
accurate across various parameter regimes.Comment: accepted for publication in Journal of Computational Physics; 61
pages, 15 figures, 11 table
Comment on Viscous Stability of Relativistic Keplerian Accretion Disks
Recently Ghosh (1998) reported a new regime of instability in Keplerian
accretion disks which is caused by relativistic effects. This instability
appears in the gas pressure dominated region when all relativistic corrections
to the disk structure equations are taken into account. We show that he uses
the stability criterion in completely wrong way leading to inappropriate
conclusions. We perform a standard stability analysis to show that no unstable
region can be found when the relativistic disk is gas pressure dominated.Comment: 9 pages, 4 figures, uses aasms4.sty, submitted for ApJ Letter
A unified hyperbolic formulation for viscous fluids and elastoplastic solids
We discuss a unified flow theory which in a single system of hyperbolic
partial differential equations (PDEs) can describe the two main branches of
continuum mechanics, fluid dynamics, and solid dynamics. The fundamental
difference from the classical continuum models, such as the Navier-Stokes for
example, is that the finite length scale of the continuum particles is not
ignored but kept in the model in order to semi-explicitly describe the essence
of any flows, that is the process of continuum particles rearrangements. To
allow the continuum particle rearrangements, we admit the deformability of
particle which is described by the distortion field. The ability of media to
flow is characterized by the strain dissipation time which is a characteristic
time necessary for a continuum particle to rearrange with one of its
neighboring particles. It is shown that the continuum particle length scale is
intimately connected with the dissipation time. The governing equations are
represented by a system of first order hyperbolic PDEs with source terms
modeling the dissipation due to particle rearrangements. Numerical examples
justifying the reliability of the proposed approach are demonstrated.Comment: 6 figure
Thermodynamic interpolation for the simulation of two-phase flow of non-ideal mixtures
This paper describes the development and application of a technique for the rapid interpolation of thermodynamic properties of mixtures for the purposes of simulating two-phase flow. The technique is based on adaptive inverse interpolation and can be applied to any Equation of State and multicomponent mixture. Following analysis of its accuracy, the method is coupled with a two-phase flow model, based on the homogeneous equilibrium mixture assumption, and applied to the simulation of flows of carbon dioxide (CO2) rich mixtures. This coupled flow model is used to simulate the experimental decompression of binary and quinternary mixtures. It is found that the predictions are in good agreement with the experimental data and that the interpolation approach provides a flexible, robust means of obtaining thermodynamic properties for use in flow models
HDG-NEFEM with Degree Adaptivity for Stokes Flows
This paper presents the first degree adaptive procedure able to directly use the geometry given by a CAD model. The technique uses a hybridisable discontinuous Galerkin discretisation combined with a NURBS-enhanced rationale, completely removing the uncertainty induced by a polynomial approximation of curved boundaries that is common within an isoparametric approach. The technique is compared against two strategies to perform degree adaptivity currently in use. This paper demonstrates, for the first time, that the most extended technique for degree adaptivity can easily lead to a non-reliable error estimator if no communication with CAD software is introduced whereas if the communication with the CAD is done, it results in a substantial computing time. The proposed technique encapsulates the CAD model in the simulation and is able to produce reliable error estimators irrespectively of the initial mesh used to start the adaptive process. Several numerical examples confirm the findings and demonstrate the superiority of the proposed technique. The paper also proposes a novel idea to test the implementation of high-order solvers where different degrees of approximation are used in different elements
Coupled, Physics-Based Modeling Reveals Earthquake Displacements are Critical to the 2018 Palu, Sulawesi Tsunami
The September 2018, Mw 7.5 Sulawesi earthquake occurring on the Palu-Koro strike-slip fault system was followed by an unexpected localized tsunami. We show that direct earthquake-induced uplift and subsidence could have sourced the observed tsunami within Palu Bay. To this end, we use a physics-based, coupled earthquake–tsunami modeling framework tightly constrained by observations. The model combines rupture dynamics, seismic wave propagation, tsunami propagation and inundation. The earthquake scenario, featuring sustained supershear rupture propagation, matches key observed earthquake characteristics, including the moment magnitude, rupture duration, fault plane solution, teleseismic waveforms and inferred horizontal ground displacements. The remote stress regime reflecting regional transtension applied in the model produces a combination of up to 6 m left-lateral slip and up to 2 m normal slip on the straight fault segment dipping 65∘ East beneath Palu Bay. The time-dependent, 3D seafloor displacements are translated into bathymetry perturbations with a mean vertical offset of 1.5 m across the submarine fault segment. This sources a tsunami with wave amplitudes and periods that match those measured at the Pantoloan wave gauge and inundation that reproduces observations from field surveys. We conclude that a source related to earthquake displacements is probable and that landsliding may not have been the primary source of the tsunami. These results have important implications for submarine strike-slip fault systems worldwide. Physics-based modeling offers rapid response specifically in tectonic settings that are currently underrepresented in operational tsunami hazard assessment
- …